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In this study, we theoretically investigate the feasibility of using laser-wakefield accelerated (LWFA) elec- 
trons for the photonuclear measurement of nuclear isomers according to the characteristics of the electrons 
obtained from LWFA experiments conducted at the Compact Laser Plasma Accelerator (CLAPA) laboratory. 
The experiments at the CLAPA show that a stable electron beam with an energy of 78-135 MeV and a charge 
of 300-600 pC can be obtained. The bremsstrahlung spectra were simulated using Geant4, which suggests that 
a bremsstrahlung source with a peak intensity of 109 photons/s can be generated. Theoretical calculations 
of isomer production cross-sections from the photonuclear reactions on six target nuclei, 197 Au, 180Hf, 159Tb, 
1151n, 103Rh, and Zr were performed and compared with the available experimental data in EXFOR, which 
suggest that further experiments are required for a series of photonuclear reaction channels. Flux-averaged cross 
sections and isomer ratios (IR) resulting from such bremsstrahlung sources are theoretically deduced. The re- 
sults suggest that IR measurements can be used to constrain nuclear components, such as y strength function 
and optical model potential. In addition, the detection of the decay characteristics was evaluated with Geant4 
simulations. The use of the LWFA electron beam and its bremsstrahlung for photonuclear studies involving 


nuclear isomers is anticipated. 
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I. INTRODUCTION 


Since the theoretical prediction by Soddy [1] in 1917 and 
the first experimental observation by Hahn [2] in 1921, nu- 
clear isomers have received intensive attention from the per- 
spective of fundamental research and applicational develop- 
ment [3, 4]. In particular, the effective population and manip- 
ulation of nuclear isomers has been an active topic in recent 
years, driven by the prospect of utilizing their intrinsic or de- 
pleted decays for the development of nuclear clocks, nuclear 
batteries, and gamma lasers [5-9]. Moreover, the exploration 
of nuclei with isomeric states in nuclear astrophysics allows 
us to determine how they affect the creation of elements in the 
universe and how they eventually contribute to the makeup of 
life in our cosmos [10]. 

In recent decades, the worldwide development of high- 
intensity and high-repetition-rate laser systems has opened 
up new research opportunities for particle acceleration [11- 
17] and nuclear photonics [18]. In particular, the ultra-short 
pulse radiation source and plasma environment generated by 
the interaction between the laser and target provide unique 
physical environments for the study of nuclear isomers. Us- 
ing a table-top hundred TW laser, Feng et al. [19] demon- 
strated that the femtosecond pumping of "Kr with a peak 
efficiency of 2.34 x 10!5 s7! can be achieved through the 
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Coulomb excitation of ions with quivering electrons during 
the interaction between laser pulses and clusters at nearly 
solid densities. Wang et al. [20] proposed a novel approach 
to populate the nuclear clock isomer 7*?"Th by combining 
laser-driven electron recollision and electronic nuclear exci- 
tation, for example, Coulomb excitation, which does not re- 
quire precise knowledge of the isomeric energy. The com- 
mon idea of these schemes is to make use of the dense low- 
energy electrons produced in laser-plasma interactions and 
the relatively high electronic excitation cross sections of ion- 
ized atoms in plasma. A comprehensive comparison of the 
possible physical processes for isomer population and ma- 
nipulation, including nuclear excitation by electron capture, 
excitation by electron transition, electron bridge, and pho- 
toactivation, was also given in Ref. [21], which indicates that 
nuclear excitations assisted by electrons would be far more 
efficient than straightforward photon excitation. This is rea- 
sonable for isomer manipulation or triggered isomer deple- 
tion circumstances, where the target nucleus must remain un- 
changed, and the excitation energy required is as low as the 
difference between the atomic electron shells. However, this 
is somehow not the case when it comes to the isomer popu- 
lation, where the isomer can be populated without extra scru- 
ples on the disintegration of the target nucleus. Specifically, 
the cross-section of photonuclear reactions can be signifi- 
cantly enhanced by giant dipole resonances (GDR) to hun- 
dreds of millibarns above the neutron separation threshold, 
whereas only the narrow discrete resonance from the isolated 
bound states is generally mentioned in most of the discus- 
sions for photo-activation of nuclear isomers. Utilizing (7,7’) 
scattering of the collective excitation mode, Pan et al. [22] 
have depicted the prospect of accumulating isomers of med- 
ical interest to a clinically applicable activity with an intense 
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Fig. 1. (Color online) Panels (a), (b), and (c) are the exemplary LWFA electron spectra obtained from the back pressure of 33, 36, and 39 
bar, respectively. (d) The bremsstrahlung energy spectrum obtained from the Geant4 simulations considering the averaged energy spectra and 
angular spread of 100 LWFA electron shots. (e) The correlation of bremsstrahlung energy against its angular spread for electrons at 33 bars. 


monochromatic y-ray beam. Meanwhile, quasimonochro- 
matic high-intensity electron pulses generated by laser wake- 
field acceleration (LWFA) in gas jets are commonly available 
in high-power laser facilities, which can be converted into 
bremsstrahlung radiation, thus providing vast research oppor- 
tunities for photonuclear studies. Although previous studies 
have utilized laser-driven photonuclear reactions for electron 
diagnosis [23], nuclear waste transmutation [24], and medi- 
cal isotope production [25], little attention has been paid to 
nuclear isomers. Consequently, a systematic investigation is 
required to evaluate the potential of the isomer population by 
laser-driven photonuclear excitation. 


In this work, we theoretically investigate the potential 
to accumulate nuclear isomers by laser-driven photonuclear 
excitation through Monte Carlo simulations based on the 
electron characteristics obtained from LWFA experiments. 
LWFA experiments at the Compact Laser Plasma Accel- 
erator (CLAPA) laboratory of Peking University are intro- 
duced, and the conversion from electrons to bremsstrahlung 
is simulated with the Geant4 [26]. Then, theoretical cal- 
culations on the photonuclear cross-sections are performed. 
A series of nuclear isomers were selected as the calcula- 
tion candidates, considering both the nuclear flow of pho- 
tonuclear reactions and the availability of stable natural tar- 
get nuclei. The corresponding photonuclear cross-sections 
for isomer generation were evaluated using TALYS [27, 28]. 
To verify the reliability of the TALYS inputs, a benchmark 
of the theoretical predictions against the available experi- 
mental data extracted from the experimental nuclear reac- 
tion database, EXFOR [29] is performed. The flux-averaged 
cross-sections and flux-averaged isomer ratios were deduced 


according to the theoretical photonuclear cross-sections and 
simulated bremsstrahlung spectra. Furthermore, the photonu- 
clear cross-sections of different reaction channels were incor- 
porated into the simulation toolkit Geant4-GENBOD [30-32] 
to estimate the feasibility of detecting the nuclear isomers of 
interest. In the simulations, the typical laser-driven electron 
beam parameters extracted from the experiments are consid- 
ered, as well as the physical processes of the radioactive de- 
cay of various reaction residuals and energy deposition in 
commercially available LaBr3 scintillation detectors. A com- 
prehensive discussion is provided in the context of the pro- 
duction and detection of the nuclear isomers of interest. 


I. LWFA ELECTRONS AND BREMSSTRAHLUNG 
CONVERSION AT CLAPA 


We conducted experiments to accelerate electrons for 
bremsstrahlung generation with the 200 TW laser facility at 
the CLAPA of Peking University, which can deliver laser 
pulses with a duration of 30 fs, energy of 4 J, and repeti- 
tion rate of 5 Hz. The laser beam was focused on a spot of 
21 x 23 um full width at maximum onto the leading edge of 
a4 x 1 mm supersonic gas jet. The effective energy on the 
target was approximately 1.8 J with a laser intensity of ap- 
proximately 1.5 x 10-!? W/cm7?. In the interaction between 
the laser pulses and the mixed gas made of 5% N2 and 95% 
He2, monochromatic energetic electrons can be produced by 
LWFA. The electrons propagate through a thin Be window; 
heading out of the vacuum chamber to the magnetic spec- 
trometer placed 100 cm downstream of the gas jet, where the 


Fig. 2. (Color online) The nuclear flow of photonuclear reactions and the mechanism for isomer generation. Here Sn denotes the neutron 


separation energy of the target nucleus 180Hf, 


deflected electron signals can be diagnosed using fluorescent 
screens and CCD cameras. The electron charge was measured 
by a Turbo Integrating Current Transformer (Turbo-ICT) as- 
sembled right out of the Be window. 


The back pressure of the gas jet was varied from 33 bar to 
39 bar, obtaining LWFA electrons with a peak energy of 78- 
135 MeV. One hundred continuous shots were recorded for 
each gas pressure to test the stability of the accelerated elec- 
trons. The repetition rate remained at 0.25 Hz to ensure that 
the vacuum level of the laser-gas interaction chamber could 
be recovered for the LWFA. Typical electron distributions 
recorded by the magnetic spectrometer are shown in Fig. 1 
(a)-(c). Electron energies of 78 + 10, 103 + 14, and 135 
+ 20 MeV were obtained at 39, 36, and 33 bar, respectively. 
The divergence of the electrons was within 4.3 mrad and the 
directional stability was better than 5.2 mrad. The electron 
charge was measured to be 300—600 pC. 


The averaged energy and divergence distributions of the 
electrons obtained at 39, 36, and 33 bar were incorporated 
into Geant4 to simulate the resulting bremsstrahlung charac- 
teristics. In the simulations, the electrons were directed to a 
2-mm-thick Ta converter 1 m away from the electron gener- 
ation point, inducing bremsstrahlung radiation. The electron 
charge was set to 300 pC, which is a conservative estimate 
for the number of LWFA electrons produced per shot. The 
bremsstrahlung photons were recorded in a plane perpendicu- 
lar to the propagation direction assembled 20 cm downstream 
from the converter. The simulated bremsstrahlung energy dis- 
tributions are shown in Fig. 1 (d). The photon intensity shows 
a typical decaying trend with the increase in photon energy, 
which can fully cover the GDR region of the selected target 
nuclei. The conversion efficiency from electrons to photons 
increases with increasing electron energy. The total photon 
flux can reach the order of 4.35 x 108 photons/shot, which 
corresponds to a peak photon flux of 1.24 x 10!° photons/s, 
considering a pulse duration spread of 35 ps inside the con- 
verter. The correlation between the photon energy and polar 
angle for the incident electrons at 33 bar is shown in Fig. 1 
(e). It is shown that the angular spread of the photons tends to 
narrow with an increase in photon energy. The angular spread 


of most of the photons is within 0.4 rad. With such a photon 
beam with a high peak flux and ultrashort pulse duration, one 
can potentially measure the nuclear isomers with short half- 
lives produced in photonuclear reactions. Accordingly, the 
reaction data covering the GDR and quasi-deuteron regions 
were probed. 


HI. PRODUCTION ROUTE AND PHOTONUCLEAR 
CROSS-SECTIONS 


A. Theory 


The nuclear flow of photonuclear reactions and mechanism 
for isomer generation are depicted in Fig. 2. When a target 
nucleus is excited above the particle separation energy the ex- 
cited compound subsequently decays through the emission of 
particles (such as neutrons, protons, and alpha particles) or 
photons to different energy levels of the residual, including 
the isomer state. At energies of 10—40 MeV the GDR result- 
ing from the out-of-phase oscillation between the proton fluid 
and neutron fluid would clearly give rise to the odds of pho- 
ton absorption and compound nucleus formation. At 40-130 
MeV, the photon interacts with the nucleus by the quasi- 
deuteron mechanism (QD), where the compound nucleus de- 
cays predominantly via the emission of several nucleons, but 
mostly neutrons. The decay of the compound nucleus, which 
would lead to the population of isomers can be well described 
by the Hauser-Feshbach reaction model, which is based on 
the fundamental Bohr hypothesis that the reaction occurs by 
means of the intermediary formation of a compound nucleus 
that can reach a state of thermodynamic equilibrium. The 
Hauser-Feshbach model averages over many resonances in 
the formation of a compound nucleus and therefore can be 
applied if the nuclear level density (NLD) in the compound 
nucleus is sufficiently high at the compound formation en- 
ergy, which is the excitation energy at which the compound 
nucleus is formed. 
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Fig. 3. (Color online) The calculated (y,xn) cross sections for 197 Au, 180Hf, 159Tp, 115 In, 103Rh, Zr. The available experimental data from 


the EXFOR database are also shown for comparison [33-53]. 
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Fig. 4. (Color online) The isomer production cross sections for the six selected target nuclei. The available experimental data from the 


EXFOR database are also shown for comparison [54—60]. 


B. Verification of photonuclear cross-sections 


Generally, photonuclear reactions tend to push the involved 
nucleus toward the down-left side of the nuclear chart. This 


places a clear limit on the possible populated isomers when 
considering the stability of the target nucleus. According 
to the nuclear flow of photonuclear reactions and the stabil- 
ity of the target nucleus, we selected a series of stable nu- 
clides, including !97Au, !8°Hf, °Tb, 1151In, Rh, and Zr 
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Fig. 5. (Color online) The FACS for (y,xn) reactions for the six selected target nuclei. The available experimental data from the EXFOR 


database are also shown for comparison [61, 62]. 


as candidates for theoretical calculations of reaction cross- 
sections. The cross sections for typical photonuclear reac- 
tions, such as (y, Y), (y, xn), (Y, p), and (y, œ), are cal- 
culated using the software package TALYS, which provides 
a complete description of all reaction channels and observ- 
ables, and many state-of-the-art nuclear models covering all 
the main reaction mechanisms encountered in light-particle- 
induced nuclear reactions are included. The nuclear structure 
components used in the TALYS calculations were as follows: 
The nuclear masses are taken from the Atomic Mass Eval- 
uation 2016 (AME2016) [67] whenever available, whereas 
the HFB-27 nuclear masses [68] are taken into account when 
the AME2016 mass data is not available. The discrete exper- 
imental levels compiled in the RIPL-3 library [69] and the 
continuum level spectrum represented by the nuclear level 
densities (NLDs) are both considered in the calculations. 
The NLDs were obtained from predictions in the Fermi gas 
model plus a constant temperature. The photon strength func- 
tions (SFs) obtained from the Brink-Axel Lorentzian model 
were used to calculate the electromagnetic transmission co- 
efficients for the photon channel. The phenomenological op- 
tical model potentials (OMPs), i. e., the local and global pa- 
rameterizations of Koning and Delaroche are employed to de- 
termine the transmission coefficient for the particle channels. 


The theoretical photonuclear reaction cross-sections based 
on TALYS for the selected targets and their corresponding ex- 
perimental data available in the EXFOR database are shown 
in Fig. 3. It can be seen that the (y, n) reaction cross- 
section is generally dominant in the GDR region, with a max- 
imum of hundreds of mb. As the excitation energy contin- 
ues to increase, multiple emission reactions take over, whose 


maximum cross-section decreases with an increase in the 
number of emitted particles. The theoretical calculations 
agree with the EXFOR data, which is obtained from pre- 
vious measurements using photon sources, such as (tagged) 
bremsstrahlung [33-39], laser-Compton y rays [40—44], and 
positron annihilation in flight [45-53]. This indicates the re- 
liability of the TALYS inputs and calculations. In addition, 
it is noticeable that the previous measurements have not fully 
covered all the photonuclear reaction channels, although the 
(y,xn) reaction data for the most studied nucleus 197 Au are 
rather abundant. Future experiments are expected to provide 
new data for these reaction channels. The isomer production 
cross-sections resulting from photonuclear reactions on the 
six selected target nuclei are shown in Fig. 4. Isomer produc- 
tion from (y, n) reactions is generally greater than that from 
the other reaction channels, with a maximum cross-section of 
hundreds of ~100 mb. It is also noted that experimental iso- 
mer production cross-sections are very rare, which can be po- 
tentially measured using state-of-the-art LCS facilities [70- 
72] in the future. 


IV. PRODUCTION YIELDS AND FLUX-AVERAGED 
ISOMER RATIOS 


A. Isomer production yields 


In bremsstrahlung-driven photonuclear measurements, iso- 
mers and their decay have long been interesting topics be- 
cause they can provide constraints for nuclear structures. Un- 
der bremsstrahlung radiation, the flux average cross sections 
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Fig. 6. (Color online) The isomer ratios of IRm2/¢ (a-c), IRm2/g (d-f), and IRm2/(m1+g) (g-i) for the 197 Au(y, n) reactions calculated with six 
sets of eight sets of SFs, two sets of OMPs, and NLDs available in TALYS 1.9. The nuclear structure ingredients used for the calculations are 
indexed by numbers and described in the main text. The experimental data from EXFOR are also shown for comparison [61, 63—66]. 


(FACS) of a given photonuclear reaction can be given by: 


or (E:)= f o(E) (EAE / f oma, w 


En En 


where o (E) is the reaction cross-section at photon energy E, 
®(E) is the photon flux, Ee is the electron energy, En is the 
photonuclear reaction threshold. 

Using the bremsstrahlung driven by laser-plasma- 
accelerated electron beams described in Sect. II, we 
conducted experimental measurements on the (y,xn) cross 
sections of !?’Au and compared them against the theoretical 
estimations based on Eq. 1. The experimental details and 
data analysis can be found in Ref. [73]. The agreement 
between experimental results and theory was found, which 
demonstrates that the procedure of measuring FACS with 
laser-driven bremsstrahlung is feasible. With this premise, 
we further evaluated the expected FACS of possible isomer 
residuals when laser-driven bremsstrahlung beams were used 
for the photoactivation of 197 Au, 0HF, 159Tb, 115Tn, 103Rh, 
and °°Zr, as shown in Fig. 5. Noticeably, for the FACS 


to produce the nuclear isomers 196ml Ay 195m Au, 179ml Hf, 
158mTh, 1l4mjn, 10lmRh and 8°™Zr were relatively larger 
than the others, exceeding the value of 10 mb. Moreover, 
the available experimental FACS to produce nuclear isomers 
for these target nuclei is very limited. The evaluated FACS 
for the production of !°6™? Au from the (y, n) reaction on 
197 Au agrees well with the experimental measurements 
made by Naik et al. [61], while those for the production 
of 102m101m,100m.9°MRh produced from (y,xn) reactions on 
'03Rh are all lower than the experimental values obtained 
by Rahman et al. [62]. The FACS for the isomers of the 
other isotopes produced from reactions, such as (y, p) 
and (y, Œœ) reactions are also shown in the Supplementary 
Material, which are less abundant in the experimental data. 
These results show that experimental efforts to measure 
FACS involving isomer residuals are expected to resolve 
discrepancies and deficiencies in photonuclear data in the 
QD region. 
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Fig. 7. (Color online) (a) The recorded y-ray spectra for the photonuclear products of °°Zr when the cooling time fe is 1 ms, 10 s, 10 m, 1 h, 
and 1 d. The data acquisition time is set to 14 days. Here the offline detection of delayed radiation is imitated by applying a time window for 
the energy deposition in the detector. (b) The dependence of recorded counts on the energy deposition and detection time. 


B. Dependence of isomer ratio on nuclear structure 
ingredients 


The ratio of the probabilities of forming the isomer and the 
ground state is called the isomeric ratio, which can be given 
by: 


omi(Es) Je? mi(E)P(E)AE / Jg? P(E)dE 


IRmi/g = = 
© o Ee) JE og(E)®(E)dE/ [E &(E)dE 


where oi (Ee) and o$, (Ee) are the FACS to produce the i-th 
isomer and ground state, respectively. Omj(E) and o,(E) are 
the production cross-sections for the i-th isomer and ground 
state, respectively. The production threshold for the ith iso- 
mer Eh differs slightly from that for the ground state Em. In 
this study, we theoretically investigated the sensitivity of IR 
to three nuclear structure ingredients, including NLD, OMP, 
and SF. In particular, the IRs of '°™Au to 1°68 Au (IRmzg); 
196m! Ay to 968 Au (Rm12), and 196m2 Au to the sum of 1°68 Au 
and 196ml Ay (IRm2/(m1+g)) produced from the (y, n) reactions 
of Au are evaluated using different sets of NLDs, OMPs, and 
SFs available in TALYS. 

The results are predicted using six sets of NLDs (1 - Con- 
stant temperature Fermi gas model [74], 2 - Back-shifted 
Fermi gas model [75], 3 - Generalised superfluid model 
[76], 4 - HFB-Skyrme model [77], 5 - HFB-Skyrme model 
with combinatorial method [78], and 6 - Temperature(T)- 
dependent HFB-Gogny model [79]), two sets of OMPs (0 
- phenomenological Wood-Saxon potential [80], 1 - semi- 
microscopic JLMB potential with HFB-Skyrme matter den- 
sity [68, 81-83]), and eight sets of PSFs (1 - Generalized 
Lorentzian [84], 2 - Brink-Axel Lorentzian [85, 86], 3 - 
HFBCS-QRPA model [87], 4 - HFB-Skyrme-QRPA model 
[88], 5 - Hybrid model [89], 6 - T-dependent HFB-Skyrme- 
QRPA model [88, 90], 7 - T-dependent RMF model [91], and 
8 - HFB-D1M-QRPA model [92]) are shown in Fig. 6. It 


is shown that the IRs grow rapidly with the increase in inci- 
dent electron energy at low energies and tend to saturate at 
great energies of tens to hundreds MeV, where discrepancies 
resulting from different models occur. Specifically, IRmi/g is 
not sensitive to the three nuclear ingredients, whereas both 
IRm2/mi+g) and IRm2/g show relatively strong sensitivity to 
the SFs and NLDs. This suggests that the existing experimen- 
tal data and further measurements aimed at IRs, especially in 
high-energy regions, can be used to constrain these nuclear- 
structure ingredients. 


V. DETECTION OF CHARACTERISTIC DECAYS 


In photonuclear experiments, there are two ways to mea- 
sure the reaction cross-sections: online detection of photo- 
induced particle emissions (neutrons, protons, and œ parti- 
cles) and offline detection of the characteristic y decays of 
radioactive residuals. The former is generally used in experi- 
ments conducted with energy-tunable, quasi-monochromatic 
probes, such as laser-Compton y-ray sources and positron an- 
nihilation in flight because the energy spectra of the emitted 
particles would exhibit easily identifiable reaction kinematics. 
In bremsstrahlung-driven photonuclear experiments, offline 
detection is typically performed. The detection of residuals 
depends not only on their decay characteristics and produc- 
tion yields but also on the energy resolution of the detector. 

We performed Geant4 simulations to evaluate the feasibil- 
ity of detecting the considered isomers produced in photonu- 
clear reactions. The simulation setup is arranged as follows. 
The energy distribution and angular spread of the bremm- 
strahlung produced by the LWFA electrons are set accord- 
ing to the experimental results obtained at a gas pressure of 
36 bar. The electron charge was set as 100 shots at 300 pC. 
The bremsstrahlung photons irradiated the activation targets 
with a thickness of 5 mm, placed 20 cm downstream of the 
converter, inducing photonuclear reactions and producing ra- 


Table 1. The expected counts, half-lives, y-ray energies, and intensities of the principal decays from the photonuclear products of !97 Au, 
180Hf, 159Tb, 1151n, !3Rh, and Zr. The cooling time fe is assumed as 1 us. The data acquisition time is set as 14 days. The level energy of 
195m? Au is labeled as 2460.84+X keV because its transition energy to the level at 2460.84 keV is not determined yet. 


Isotope Half-life Level energy (keV) Decay energy (keV) Intensity (%) Net counts (14 d) 


DMA 7.73s 409.16 279.01 70.9 5.94 x 107 
19% Au 6.1669 d 0.00 355.73 87 4.44 x 10° 
196ml Au 8.15 84.66 9.71 25 N/A 
196m2 Au 96h 595.66 147.81 43.5 2.02 x 102 
5 au 186.01 d 0.00 98.86 11.21 4.03 x 102 
195ml Ay 30.55 318.60 261.75 68.7 1.69 x 10° 
195m2 Au 12.89 us 2460.84+X X N/A N/A 
19% Au  38.02h 0.00 328.46 60.4 1.23 x 10° 
194ml Au 600 ms 107.40 45.32 5.4 N/A 
194m2 Au 420 ms 475.80 170.78 27.3 2.22 x 102 
193 Au  17.65h 0.00 255.57 6.5 3.52 x 103 
I93mAy 39s 209.203 257.97 67.14 2.22 x 104 
180mHf  5,53h 1141.50 332.27 94.00 N/A 
179m Hf 18.67 s 375.00 214.34 95.3 1.29 x 106 
179m2Hf 25.05 d 1105.70 453.59 68 2.77 x 10° 
l8mlyf 4s 1147.40 426.36 97 1.73 x 104 
M8m2Hf 31y 2446.10 574.22 88.6 7.25 x 102 
I7ml yf 1.09s 1315.50 208.37 73.1 1.24 x 103 
Wm2Hf 51.4m 2740.00 277.30 75 4.42 x 102 
158Tb 180 y 0.00 944.19 44.40 5.65 x 102 
158mi Th 10.75 110.30 110.30 0.92 1.37 x 10° 
158m2Th 0.40 ms 388.40 171.07 50.4 2.76 x 103 
S7Th  7ly 0.00 54.50 0.0084 N/A 
156Tb 5.35 d 0.00 199.19 67 6.35 x 10° 
omlTh 24.4h 49.60 49.63 74.1 N/A 
156m2Th 5.3h 88.40 88.40 1.15 8.59 x 102 
ISTh  5.32d 0.00 105.32 25.1 4.02 x 102 
1Smqy 4.486 h 336.20 336.24 45.90 1.59 x 10° 
Mw 719s 0.00 1299.83 0.14 2.43 x 103 
M4mlf, 49.51 d 190.30 725.24 4.4 2.66 x 104 
14min 43.1 ms 501.95 311.65 89.85 1.05 x 10° 
113mJn 99.476 m 391.70 391.70 64.94 1.68 x 10° 
2%) 14.88 m 0.00 617.52 1 9.92 x 103 
112mlīn 20.67 m 156.60 156.61 13.33 1.15 x 104 
112m2%y 2.81 us 613.82 187.93 N/A 1.13 x 104 
liin 2.8047 h 0.00 245.35 94.1 7.20 x 104 
limn 77m 537.00 537.00 87 9.99 x 103 
1035mRh 56.114m 39.80 39.76 0.07 N/A 
12Rh 207.3d 0.00 475.06 46 1.57 x 10° 
102mRh 3.742 y 140.70 697.49 44 1.18 x 10° 
IRh 33y 0.00 198.01 73 5.60 x 103 
10ImRh 434d 157.30 306.86 81 9.21 x 10° 
100Rh =. 20.8h 0.00 539.51 80.6 2.91 x 10° 
100mRh 4.6m 107.60 539.60 1.68 2.91 x 10° 
"Rh ě 161d 0.00 528.24 37.9 1.58 x 104 
99mMRh 47h 64.60 340.80 72 4.80 x 10° 
90m7, 809.2 ms 2319.00 2318.96 83.00 1.28 x 10° 
897, = 78.41h 0.00 909.15 99.04 3.37 x 10° 
89m7; 4.161 m 587.80 1507.40 6.06 1.43 x 105 
88Zr 83.4d 0.00 392.87 97.29 5.20 x 104 
877,  1.68h 0.00 1227.00 2.8 5.03 x 102 
87m7Zr 14s 335.80 201.02 96.4 6.32 x 103 


867,  16.5h 0.00 242.80 95.84 1.38 x 104 


dioactive residuals. The y decays of the radioactive residu- 
als, including the isomers, were recorded by a LaBr3 detector 
with dimensions of ọ 8 cm x 10 cm assembled at 90° with 
respect to the beam propagation. The energy resolution is set 
as 3%, which is achievable for commercially available LaBr3 
detectors. To simulate isomer production and decay radia- 
tion detection in a one-stop manner, we adopted the Geant4- 
GENBOD [30-32] toolkit and further implemented the ra- 
dioactive decay process of the radioactive residuals (includ- 
ing the isomers). The photonuclear reaction cross-sections of 
the related isotopes were calculated with a routine similar to 
Sect. II B. 


An exemplary correlation of energy deposition versus de- 
tection time for the radioactive residuals induced by photonu- 
clear reactions on Zr is given in Fig. 7. In activation exper- 
iments based on laser-plasma acceleration, there is a period ( 
cooling time t.) for opening the reaction chamber before the 
activated samples are handed to the detector for data acqui- 
sition. If fo is larger than 1 h and 10 s, the decay informa- 
tion of the products ®°™Zr (T, /2 = 4.16 m) and iu a (T) /2 
= 809 ms) will be lost. To detect isomers with short half- 
lives, such as °°™Zr the detector is placed near the activation 
target in an "online" manner. The backgrounds induced by 
the backscattered bremsstrahlung, electrons, and electromag- 
netic pulses can be cut off by applying a time window with a 
delay signal generator. The correlation of energy deposition 
versus detection time for the radioactive residuals induced by 
the photonuclear reactions on all six selected targets is pre- 
sented in the Supplementary Material. The principal decay 
characteristics and expected counts within three half-lives for 
the isomers produced from the photonuclear reactions of the 
six selected isotopes are shown in Table 1. The photonuclear 
products 196ml Au, 194ml Au, 157Tb, 156ml Th, and 103mph are 
not measurable for the LaBr3 detector because their decay 
energy is too little and the attenuation in the target is strong. 
The measurement of !°™?Au is unlikely because its decay 
energy and intensity has not yet been determined. For !*°Hf, 
the decay of '8°™Hf is difficult to measure because the (y, y) 
reaction cross-section is very low. In addition, most of the 
isomer residuals can be measured with a net count of 10°~° 
within 14 days of data acquisition. The decay energies of 
some products, such as 1° Au and !93™Au are very close and 
cannot be disentangled from the recorded energy spectrum of 
a detector with poor energy resolution. However, this can be 
remedied by applying a fit to the recorded counts in a few 
data-acquisition cycles. 
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